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The construction is carried out explicitly in 2-d, on an arbitrary orientable 
manifold without boundary. It can be easily converted into a computer code. 
The equivalence, on a sphere, of Majorana fermions and Ising spins in 2-d is 
rederived. The method can, in principle, be extended to higher dimensions. 
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1 Introduction 



1.1 Preamble 

The statistical mechanics of random manifolds is being intensely studied since 
a dozen years (see the reviews ) . The interest of endowing these manifolds 
with fermionic degrees of freedom seems rather evident. It turns out that 
little has been done in this direction. A notable exception is the paper by 
Bershadsky and Migdal @, where it is demonstrated that the Ising model on 
a fixed planar 2-d graph is equivalent to a certain Majorana fermion theory 
on the dual graph. However, we do not find the discussion of ref. [|2| fully 
satisfactory. The authors ignore the covariance aspects of the problem and 
therefore several features of the model, which have a natural explanation, 
come out as "miracles". Furthermore, it does not seem possible to extend 
their arguments beyond 2-d. Some useful ideas can be found in the pioneering 
papers of the Columbia group, in particular in ref. M , but they do not discuss 
at all the topological aspects of defining a spin structure on a piecewise linear 
manifold. 

The aim of the present paper is to propose a rather general method of 
constructing the Wilson fermion action on a piecewise flat manifold, made 
up by gluing equilateral simplices at random. The metric is assumed to 
have the Euclidean signature. Our motivation for studying this problem 
originates from our involvement in the study of simplicial gravity. In this 
context, limiting oneself to such manifolds is a common practice and seems 
justified by the results obtained in the study of non-critical strings: for a 
class of exactly solvable models in 2-d it has been shown that the continuum 
and the discrete version belong to the same universality class, when the 
"dynamical triangulation" recipe is adopted |T[]. According to this recipe 
the sum over metrics in the Feynman integral is indeed replaced by the sum 
over all triangulations of the above mentioned type. In order to make the 
discussion as clear as possible we shall focus on 2-d. The results of ref. 
will, of course, be recovered. In the last section we shall briefly discuss the 
extension of the idea to arbitrary triangulations and to higher dimensions. 

The reader will notice that the concept of a continuous manifold appears 
in our discussion as a scaffolding that helps to achieve the goal of putting a 
spinor field on an abstract graph, with specific intrinsic symmetries. 
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1.2 General strategy 

How to put a spinor field on a curved manifold is explained in textbooks of 
general relativity (see e.g. Q). The key concepts are local frames, parallel 
transport and spin connection. The application of the general recipe to the 
case of a manifold discretized a la Regge is, however, not quite trivial, since 
one cannot limit the discussion to infinitesimal displacements and rotations. 

For pedagogical reasons let us recall some differential geometry. With 
each point x of the manifold is associated an arbitrary orthonormal local 
frame defined by a set of d vectors (x) called vielbeins. Here, the latin 
index j refers to the axes of the local frame transforming under local rotations 
belonging to SO(d). Using vielbeins one can eliminate, by contraction, the 
space (here Greek) indices, in order to deal only with objects transforming 
under SO(d), like in flat space. One can then introduce spinor fields, because 
the group SO(d) has spinor representations, contrary to the group GL(d) 
induced by general coordinate transformations. The program is incomplete, 
however, until one tells how to compare local frames at distinct points of 
space. When e?Jx) is parallel transported from x to x + 5x one gets a 
frame rotated with respect to the one that has been chosen at x + 8x. This 
observation is formally expressed by the equation 

efrx + Sx) = efrx) + r^(x)e{(x)5x» - <(x)ej(z)^ . (1) 

Here, is the Christoffel symbol and u) 3 ku is the spin connection, which 
takes care of the relative rotation of the neighboring local frames. Eq. ([I]) 
can be used to define a covariant derivative, which gives zero acting on the 
vielbein field (and therefore also on the metric g^ v = elejv)- 

Since we are interested in piecewise flat manifolds, we can assume that 
there exists an orthonormal reference frame in a region comprising the points 
x and x + 8x. Choosing this frame one gets from (|l]) 

e^x + 5x) = [5{ + uUx)5x"] e^x) . (2) 

The operation on the right-hand side is an infinitesimal rotation. The space 
index \x is inert since, by assumption, the parallel transport from x to x + 5x 
is trivial. 

Consider now a closed curve C in our piecewise flat manifold. It is as- 
sumed that C can be covered with a set of flat patches P a ,a = l,...,n: 
P a (lP a+ i is non-empty and there exists an orthonormal reference frame com- 
mon to the whole of P a U P a +i- The conventions are such that P n+ i = P\. 
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Figure 1: The set of patches covering a closed curve (not shown explicitly). The region 
within two overlapping neighboring patches is flat. 

We associate a common e^(a) with all points of P a . The analogue of (@) is 

4(b) = Ul(ba)el(a), (3) 

where P\, is the patch next to P a and U(ba) is a finite rotation from a to b. 
In the following, the space indices, like \i above, which are irrelevant to our 
discussion will not be exhibited. It will be convenient to eliminate the latin 
indices too, using the matrix notation to write, for example 

e(6) = U(ba)e(a) , (4) 

instead of (0). As already mentioned, the choice of the local frames e(a) and 
e(b) is arbitrary. Thus U(ba) is defined up to a local gauge transformation 

U{ba) -> G(b)U(ba)G-\a) . (5) 

A vector v(a) is uniquely defined in P a by its components along the d 
vielbeins e(a). Eq. (Q) tells how these components change as one goes from 
a given patch to its neighbor. 
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For any two neighboring patches the rotation depends on the gauge choice 
at these patches. It is easy to see that for the closed path P 1 P a — > 

P a+ i... — > Pi the global rotation matrix U(C) has the gauge transformation 

U(C) ^ G(l)U(C)G-\l) (6) 

and therefore |Tr[/(C) is a gauge independent geometrical object, whose 
deviation from unity is a measure of the curvature of the manifold. An 
orthonormal frame common to all patches does not, in general, exist and the 
chain of spin connections has inherited information about the curved metric. 

The transformation matrix in (^) belongs to the vector representation of 
SO(d). A similar transformation law holds for a spinor. The only differ- 
ence is that the corresponding transformation matrix belongs to the spinor 
representation of SO(d). 

In this paper, the simplices of the triangulation play the role of the flat 
patches above. Two neighboring simplices have a common (d— l)-dimensional 
face. The spinor field lives on simplices or, stated differently, on the vertices 
of the dual lattice. Only spinors belonging to neighboring simplices are di- 
rectly coupled. There always exists an orthonormal frame common to two 
neighboring simplices. Therefore the coupling is given by the corresponding 
spinor connection, which can be constructed without too much effort for a 
pair of simplices. The problem is that one has to define the spinor connection 
consistently all over the lattice, which is non-trivial (there exist topologies 
where this is impossible). The whole construction will be done explicitly, in 
2-d, in the next section. 

2 Two-dimensional manifolds 

2.1 Spin connection in the vector representation 

We consider a triangulation of an orientable 2-d manifold. By convention, all 
triangles are oriented counterclockwise. In particular, given a triangle a, the 
angle between e 1 (a) and e 2 (a) is tc/2 when measured counterclockwise. We 
focus on two triangles, say a and b, sharing a common link. It is convenient 
to introduce two auxiliary local frames, f(ab) and f(ba), attached to this link 
and rotated by the angle tt with respect to each other: f x {ab) (resp. f 1 (ba)) 
is perpendicular to the link and points towards the exterior of a (resp. b). 
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Figure 2: The local frames e(a) and e(b) in two neighboring triangles a and b are shown. 
The two auxiliary frames f{ab) and fiba) are also exhibited. They all have the same 
orientation. The oriented arc in the triangle a represents the angle 4> a —>b- Rotation by 
this angle brings the frame e(a) into f{ab). Similarly, the arc in b represents the angle 
(j>b-, a - The rotation by this angle brings e(b) into f(ba). Notice, that the angle is always 
measured from the local to the auxiliary frame. This is why an inverse rotation appears 
in eq. (|). 

The operation performed by the spin connection matrix U(ba) can be defined 
by the following chain of rotations 

U(ba) : e(a) -> f(ab) -> f(ba) -> e(6) . (7) 

Denote by R{4>) the rotation by angle 0: R{4>) = exp (e0) , where e is the 
rotation generator represented by the antisymmetric matrix e±2 = — e 2 i = 1. 
Hence 

U{ba) = J R- 1 (0^ a ) J R(vr) J R(0 a ^ 6 ) . (8) 

Here 4> a ^b is the angle between e 1 (a) and / 1 (a6), while (fib^a is the angle 
between e l (b) and / 1 (&a) (see Fig. 2). It should be kept in mind that these 
angles are oriented. Notice, that 

U(ba)U(ab) = 1 . (9) 

Repeating the above argument one eventually associates a connection 
matrix with every oriented link of the dual lattice. Let L n be an elementary 
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loop of the dual lattice going through n triangles labeled 1,2, ...n and let 
U(L n ) be the parallel transporter around L n . One has 



n 



Tr U(L n ) = Tr J] R(*)R{<f>. 



(10) 



k=l 



using the cyclic labeling convention (0 = n and n+1 = 1). Each factor above 
corresponds to a rotation by the angle n + 4>k-,k-i — 0fc-»fc+i equal, modulo 
27r, to 7r/3 (resp. — 7r/3) for L n oriented clockwise (resp. counterclockwise). 
Hence 



independently of the orientation of the loop. For n — 6, when the angular 
defect is zero and therefore the lattice is locally flat, the right-hand side 
equals unity, as expected. 

2.2 The spinor case 

The logic underlying the construction of the spinor connection, transforming 
a spinor field from one local frame to another, is close to that of the preceding 
subsection. There are, however, extra complications due to the fact that 
spinor rotations by <fi and (f) + In are not equivalent. The resulting sign 
ambiguities require care. 

The analogue of U(ba), to be denoted V(ba), is a matrix belonging to the 
spin 1/2 representation of the rotation group. In 2-d there are two hermitian 
Dirac matrices 7 1 and 7 2 , satisfying the usual anticommutation relations 
{7^, 7 fc } = 28^ . Another standard matrix is 75 = — «7 1 7 2 , which in 2-d is 
equal to the generator of rotations [j 1 , r y 2 ]/2i. 

The matrix representing a rotation by an angle is 



The sign ambiguity cannot be resolved unless one is in a position to control 
the angle <fi in the range to Air. The connection matrix V(ba) has the same 
sign ambiguity. In analogy to (Q) write : 



and gba is a sign factor. The parallel transport of the spinor from a to b and 




(11) 



S(0) = ±exp(* 75 0/2). 



(12) 
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Figure 3: Six possible ways a dual lattice loop can go through a triangle. The line 
segment pointing up from the centre of the triangle is the vertex flag indicating the gauge 
choice, i.e. the direction of the vielbein e 1 (k). The dashed line is the loop slightly displaced 
to the right. The sign factor z^-i k k+i is negative when the dashed line crosses the flag. 

back does not introduce any change : 

V(ab)V(ba) = 1 (14) 

which implies that 

9ab = ~9ba ■ (15) 

Replacing R by S on the right-hand side of eq. ( |T0D and including the 
sign factors one finds the spinor analogue of eq. (fTTf ) : 

n 

Tr V(L n ) = Tr J] (/ fcfc+1 5(7r)5 , (^ fc _ 1 )5 , - 1 (0 fc _ +fc+a ) . (16) 
k=l 

As in the vector case, the three rotations following the n symbol cor- 
respond to a single rotation by ±7r/3 modulo 2tt. But, since we are now 
working in the spinor representation, the 2ir is not innocent since it yields 
an extra negative sign. Hence, in general, the rotation is ± exp(±i7 5 7r/6). 
The sign in front of the exponential has to be determined carefully. It does 
not only depend on the way the dual lattice loop goes through the triangle 
k but also on the choice of the gauge, i.e. on the direction of the vielbein 
e 1 (k) (in 2-d specifying the direction of a single vielbein suffices to fix the 
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Table 1: The sign and the Kac-Ward factors. 



Fig 3 




4>k^k+l 




a 


tt/3 


5tt/3 


+ exp(-27 5 7r/6) 


b 


7T 


tt/3 


— exp(— 27 5 7r/6) 


c 


5tt/3 


7T 


— exp(— 27 5 7r/6) 


d 


5tt/3 


tt/3 


- exp(+i7 5 7r/6) 


e 


tt/3 


7T 


+ exp(+i7 5 7r/6) 


f 


7T 


5tt/3 


+ exp(+z7 5 7r/6) 



local frame). Since Tr V(L n ) is gauge invariant, we can fix the gauge at our 
convenience. We shall assume that in each triangle e 1 points from the center 
of the triangle towards one of the vertices^. Fig. 3 illustrates the six possible 
ways the loop can take through the triangle k. The direction of e l (k) is also 
indicated, and appears as a flag associated with the vertex of the dual lattice. 
The result of the calculation of the rotation matrix is given in Table 1. 

In the first three cases the path in the figure turns left. The corresponding 
elementary loop goes counterclockwise. In the remaining cases it goes clock- 
wise. The exponential factors are known as the Kac-Ward factors. Drawing 
a dashed line parallel to the loop, on the right of it, as shown in Fig. 3, one 
can see from the Table 1 that the sign factor is negative when the dashed 
line crosses the flag. Otherwise it is positive. Denote the sign factors by Zkji, 
where the letters k,j, I refer to the three successive triangles on the loop. 
Collecting all the sign factors we finally get 

^Tr V(L n ) = F(L n ) cos(mr/6) , (17) 

where 

F(L n ) = g 12 Zl23 923 z 234 ■ ■ ■ 9nl z n\2 ■ (18) 

Notice, that the loop sign factor can be defined for any closed loop, not only 
for an elementary one. 

2 Actually, the result of the calculation would be the same if the vielbein were rotated 
forth or back by an angle less than tt/3. 
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For n = 6 the sign factors on the right-hand side of flTB|) should combine 
to give F(Lq) = — 1 , because the parallel transport is trivial in flat space. 
Eq. (|T7]) is a relation between an invariant measure of the local curvature, 
on the left-hand side, and the angle deficit 2ir —mr/3, which determines the 
argument of the cosine function on the right-hand side. It is clear, that if 
one could change the deficit angle continuously, like e.g. at the top of a 
cone, the factor in front of the cosine would not change discontinuously. On 
a dynamically triangulated surface the quantization of the deficit angle is a 
lattice artifact, devoid of any deep physical significance. Therefore we set 

F{L n ) = -1 (19) 

for all elementary loops L n , whatever n is. We shall prove now that these 
constraints can be satisfied on an orientable 2-d manifold. 



2.3 Satisfying the constraints on an orientable 2-d ma- 
nifold 

We have already introduced the vertex flags. Now, we also attach flags to 
links: we put a flag on the right hand side of the dual lattice link going from 
a to b if g a b = —1. There is nothing fundamental in these flags. They are 
merely a convenient tool helping to set the spin structure on the simplicial 
manifold. 

Let us follow a counterclockwise oriented elementary loop staying always 
slightly outside of it. Since each crossed flag corresponds to a negative sign 
factor, it follows from the rule established in the preceding section that eq. 



(|T9| ) is satisfied provided we cross an odd number of flags. Following a clock- 
wise oriented elementary loop staying always slightly inside of it, one also 
requires that the number of crossed flags is odd. The two requirements are 
equivalent, since the total number of flags attached to a loop is even. In 
short, there should be an odd number of flags both outside and inside every 
elementary loop independently of its orientation. A portion of a dual lattice 
with flags put on vertices and links is shown in Fig. 4. It should be clear 
at this point of the discussion that the exact direction of a flag is irrelevant. 
One can rotate them as long as one does not cross a link of the dual lattice. 
Let us briefly outline the strategy adopted to prove that the constraint 



( jT9| ) can be satisfied on every orientable manifold: First, by direct inspec- 



tion, we check that it can be satisfied on a minimal sphere. Then, we show 
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Figure 4: A portion of a dual lattice with flags put on. The number of flags inside each 
elementary loop is odd. Notice that the dashed line following a non-elementary loop also 
crosses an odd number of flags. 



that it can be preserved when one deforms locally the geometry with an er- 
godic move, whose repeated application enables one to construct an arbitrary 
sphere. Since the constraint is satisfied for the initial configuration and it is 
preserved by the moves, it can be satisfied on an arbitrary spherical lattice. 
Higher genus surfaces can be produced by gluing spheres. We show that this 
gluing can also be done preserving the constraint. 

It is easy to convince oneself that one can satisfy the constraint ( |T§D on a 
minimal sphere, the one made up of four triangles (see Fig. 5a). A sphere of 
arbitrary size can be constructed by using the moves introduced, for example, 
in ref. [[J. One move consists of splitting an elementary loop of the dual 
lattice by inserting a new link. The inverse move consists in removing a link. 
We show below that the construction of a sphere of arbitrary size can be 
done preserving (|l^) for every elementary loop. 

When a link is added, the flags preexisting inside the loop are partitioned 
among the two new loops. Since their total number is odd, one loop gets 
an even number of old flags and the other an even number. One has to put 
five new flags, three on the new links and two on the new vertices. This can 
always be done so as to have at the end an odd number of flags inside the 
two new loops, without modifying the outside. A similar argument holds for 
the inverse move. One checks first that the number of flags to be removed 
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Figure 5: (a) A possible assignement of flags on a minimal sphere, (b) The shaded 
areas represent two identical spheres with one loop cut out. We show a way of gluing 
them together with a bridge made up of new links. Only the new flags on the bridge are 
exhibited. In each new loop there is an even number of flags, old and new, which are not 
shown explicitly. Repeating the operation enables one to create a sphere with handles. 

inside (or outside) the new loop to be created is odd (it cannot be 0!). If this 
is not the case one flips one of them outside (inside), flipping simultaneously 
inside (outside) a neighboring link flag, the one not to be removed. The new 
flag assignment satisfies the constraint if the old one did. The total number 
of flags inside two adjacent elementary loops is even. Removing a link one 
erases an odd number of flags. The final number of flags inside the new, 
larger loop is therefore odd, as it should be. 

One can extend this result to a sphere with handles. One starts by cre- 
ating the required topology and then one uses the split-join moves and the 
ergodicity argument, as above. First, duplicate a sphere and choose a pair of 
identical loops, one on each sphere. Join two such twin loops with a minimal 
number of new links as in Fig 5b. By symmetry, there is always an even 
number of old flags in the newly created loops. It is easy to convince oneself 
that new flags, three for each new link, can always be put so as to have an 
odd number of flags in the new loops. A surface with an arbitrary number 
of handles can be obtained by repeating this operation. Of course, spheres 
have to be glued in more than one place. 

Finally, let us observe that the satisfaction of the constraint (plj) by all 
elementary loops implies that it is also satisfied by an arbitrary contractable 
and self-avoiding loop C. This results from the fact that such a loop can be 
obtained by fusing elementary loops with the same orientation. The total 
number of flags inside two loops that fuse is even and one drops an odd 
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Figure 6: Fusion of two loops. The number of flags on the common line is 
always odd. 

number of flags (see Fig. 6). Hence 

F{C) = -1. (20) 

2.4 Dirac-Wilson operator, fermion loops and the Ising 
model 

The Dirac operator is defined by contracting the connection matrix V(ab) 
with the Dirac matrix 

lab = f\ab) ■ 7 = S-\<P a ^ b ) 7 1 S(<P a ^ b ) (21) 

Using eq. fllTf ) one finds that it has the following form : 

D(ab) = lab V{ab) = g ab S" 1 ^) 7 1 S(n) S(<j) b ^ a ) . (22) 

From here on we shall rather use the Dirac-Wilson operator for which we 
keep the same symbol : 

D(ab) = g ab S-\<p a ^ b ) P S(tt) S(<f> b ^ a ) , (23) 

where P = (1 + 7 : )/2 is a projector. We recall that g ab has been set following 
the rules discussed in the sect. 2.2. 

The Dirac-Wilson operator fulfills the following two conditions : 

D(ab) D(ba) = (24) 
13 



and 

D T {ab) = C D{ba) C' 1 . (25) 

where C is the charge conjugation matrix^. The former one means that the 
Wilson fermions do not propagate back and forth on the same link. The 



latter one results from eq. (|T^). For Majorana fermions one obtains from 

^ t/j{a) D{ab) 4){b) = 4>{b) D{ba) V>(a) , (26) 

which means that fermion lines are not oriented. 

Let us calculate the contribution to the partition function of a closed loop 
C = {l,2,...,n,l} : 

^(1) D(l, 2) ^(2) m D(2, 3) ^(3) . . . ^(n) D(n, 1) ^(1) 
= -Tr D(l,2) D(2,3)... D(n, 1) 

(27) 

= -F(C)Tr n- = i^ 
= -F(C)( v / 3/2) n . 

where Kj = exp (±27571/6) is the Kac-Ward factor at the j-th turn of the 
loop. The presence of the projector P makes the final result independent of 
the ± sign in the exponent of Kj. We see that the right-hand side of (|27D is 
positive when F(C) = — 1. We have already proved that this is the case for 
contractable, self-avoiding loops. Now, the loops are necessarily self-avoiding 
for Majorana fermions. From eq. fl27D easily follows the isomorphism, on a 
triangulated spherical surface, of the Majorana fermion model and of the 
Ising spin model in 2-d. It rests on the identity of the pattern of phase 
boundaries, in the Ising model, and that of closed fermion loops. We refer 
the reader to ref. for more details. 

In two dimensions and for Majorana fermions eq. (|23D can be rewritten 
in another, particularly elegant and suggestive form. Choose the represen- 
tation : 7 1 = er 3 and j 2 = a±, where cr, denotes the Pauli matrices. Then 
75 — 02 and C = 102 = e. Thus for Majorana fermionsQ ip = e ip T and 

3 C~ VC = ~¥ T , CtC = 1 and C T = -C. For Majorana fermions ip = C$ T . 
4 Using indices ip — > tp a , ip —> ip a the two equations can be rewritten in the standard 
form ipP = ip a e a P and ip a = e a pipP . 
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ip = ip e. With these conventions it is easy to check that fl23|) can be written 

D(ab) = g ab s((f) a ^ b ) ® s(0 6 ^ a ) , (28) 

where 

(cos | \ 
s = s T e = (-sinf cos |) (29) 

With our choice of the gauge the angles <p a (n), n = 1, 2, 3 are tc/3, n and 5n/3 
(the links emerging from a are ordered counterclockwise, starting with the 
flag). Hence, with each dual link n emerging from a one can associate a 
"spinbein" s n : 

'-(!) a ^(~T) « 3o » 

One has 

= s 2 s 2 = s 3 s 3 = , . 

s^ 2 = s's 3 = s 2 s 3 = y/3/2. { } 

Eq. (p8|) can be rewritten as 

£>(a&) = g ab s n(a) ® s m(b) , (32) 

where n(a) and m(6) refer to the link ab, but labeled according to the gauge 
chosen at a and b respectively. Finally, we write 

^(o) D(ab) 4,{b) = g ab $(a) s n(a) ] [^^(6)] (33) 



It should be clear that the particular choice ( j30|) is not important, only the 
invariant relations ([31]) matter. Notice, that in the above formulation one 
only works with the dual lattice, decorated with flags, with spinbeins living on 
the links of the graph. One is no longer referring explicitly to the underlying 
continuous manifold. 



2.5 A sketch of the computer implementation 

The fermion action can be constructed explicitly for use in computer sim- 
ulations. The crucial part of the construction is to choose the gauge and 
determine the dual link sign factors consistently. Hence, on each particular 
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random lattice one first puts the flags on links and vertices. In practice, 
this is most simply done recursively: one starts with an arbitrary dual loop, 
one puts an odd number of flags in its interior, one enlarges the domain with 
flags put on by considering a neighboring loop and so on. Our general results 
insure that one has no problem arriving to the last loop. 

With each dual link one associates one of the nine matrices defined by 
(|23|) , or what amounts to the same, by (p8|). These matrices can be calculated 
beforehand and stored in computer's memory. Which matrix is associated 
with a given link depends on the gauge choice at the ends of the link (i.e. 
the orientation of the vertex flags). For example, take the spinbein formal- 
ism: the flag at a given dual vertex determines the labeling of spinbeins at 
this vertex. Since each link connects two vertices, there are two spinbeins 
associated with each link and the right matrix is the one which is equal to 
their Cartesian product, see eq. fl3"2"|). 

Let L denote the number of dual links. The Dirac- Wilson operator for 
the full lattice is an L x L matrix made up of 2 x 2 matrices. The off-diagonal 
ones have the form (p2|) and are found as explained above. The diagonal ones 
are unit 2x2 matrices multiplied by a mass coming from the mass term. 



3 Conclusion 

Since the conceptual difficulty in defining a spin structure on a random man- 
ifold is associated with the sign ambiguities, the extension of our discussion 
to an arbitrary triangulation is, in principle, straightforward. When the tri- 
angles making up the lattice are not equilateral, the angles in the Kac-Ward 
factors are no longer equal to ±vr/6, but vary along the dual loop. On the 
other hand, the sign factor is a topological property of the loop and is inde- 
pendent of the triangulation details. The Dirac- Wilson operator can still be 
defined by eq. (^) but its computer implementation becomes more tedious 
because of the additional freedom in the angles entering the spin connection. 
With the definition (£23|) the factor (y / 3/2) n in fl2~T|) is replaced by the product, 



along the loop, of the cosines of the Kac-Ward angles. This no longer corre- 
sponds to a simple Ising model with a constant coupling between neighbors. 
The isomorphism between the fermion theory and the simple Ising model 
is recovered, however, if following ref. [[§ one introduces in the definition 



( 23|) of the Dirac- Wilson operator an appropriate link- depending factor. In 



substance, this amounts to define the Dirac operator by eq. (|32|) . 
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Another issue is the spectrum of the Dirac- Wilson operator. This very 
interesting problem, completely beyond the scope of this paper, deserves 
a separate study. As is known from the accumulated experience with the 
lattice fermions, spurious zero modes do occur. This pathology is likely to 
be also present for fermions interacting with the geometry. An extension of 
our discussion incorporating the recent progress in the lattice fermion theory 
would, of course, be welcome. 

We have illustrated our method of putting Wilson fermions on a randomly 
triangulated manifold by carrying the program in detail in the 2-d case. We 
do not see any major conceptual problem in extending this construction to 
higher dimensions, although we have not done this explicitly up to the end. 

Let us briefly sketch how this could perhaps be done in 3-d. The aux- 
iliary frames, f a b and fb a , have their 1st axis perpendicular to the common 
face (triangle) of the tetrahedra a and b. One can assume, without loss of 
generality that the 2nd axis is common to these two frames, which are then 
related by a rotation by ir around the 2nd axis. This does not specify these 
frames yet, because the orientation of the 2nd axis can be arbitrary. The 
natural choice is to assume that the 2nd axis points from the center towards 
one of the vertices of the triangle. There remains the freedom to make rota- 
tions by 2?t/3 within the triangle. A convenient choice of the gauge consists 
in associating the local frame e with one of the faces, in much the same way 
as for the auxiliary frames, in analogy to what was done in 2-d. 

Since the rotation group is now non-Abelian, manipulating rotations is 
less straightforward than in 2-d. However, within a tetrahedron the rotations 
relating the various frames, we have just mentioned, satisfy a simple algebra. 

The construction of the spin connection follows the rules formulated in the 
preceding section. The real problem, as in 2-d, is the consistent determination 
of the sign factors. The sign associated with a particular way an elementary 
dual loop goes through a tetrahedron is related to the fact that in the spinor 
representation the rotation connecting two faces directly can differ by a sign 
from the compound rotation, where one goes from the face to the local frame 
and then from the local frame to the other face. The bookkeeping is more 
complicated than in the 2-d case. 

The bookkeeping starts being really complicated as one attempts to fix 
the signs all over the lattice. The right strategy seems to be that already used 
in 2-d but the explicit implementation seems much more tedious. But these 
seem to be merely technical problems. We believe that the main conceptual 
issues are well illustrated by the present paper. 
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